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CN Abstract 

' — I We report a search for new gravitational physics phenomena based on Riemann- 

Cartan theory of General Relativity including spacetime torsion. Starting from the 
parametrized torsion framework of Mao, Tcgmark, Guth and Cabi, we analyze the 
!• motion of test bodies in the presence of torsion, and in particular we compute the cor- 

. !^ rections to the perihelion advance and to the orbital geodetic precession of a satellite. 

We consider the motion of a test body in a spherically symmetric field, and the motion 
H of a satellite in the gravitational field of the Sun and the Earth. We describe the torsion 

field by means of three parameters, and we make use of the autoparallel trajectories, 
which in general differ from geodesies when torsion is present. We derive the specific 
approximate expression of the corresponding system of ordinary differential equations, 
which are then solved with methods of Celestial Mechanics. We calculate the secular 
variations of the longitudes of the node and of the pericentcr of the satellite. The com- 
puted secular variations show how the corrections to the perihelion advance and to the 
orbital de Sitter effect depend on the torsion parameters. All computations are per- 
formed under the assumptions of weak field and slow motion. To test our predictions, we 
use the measurements of the Moon's geodetic precession from lunar laser ranging data, 
and the measurements of Mercury's perihelion advance from planetary radar ranging 
data. These measurements are then used to constrain suitable linear combinations of 
the torsion parameters. 
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1 Introduction 

Einstein's theory of General Relativity (GR) successfully describes gravitational physics in 
the solar system. Its predictions have passed a wide variey of precision experimental tests 
carried out in weak-field and slow-motion regime with natural bodies and artificial satellites 
[l] , p] . These tests include the measurement and verification of (quoted in terms of relative 
experimental uncertainty): Mercury's perihelion advance at 10~^ level through planetary 
radar ranging the redshift of spectral lines of a hydrogen-maser frequency standard 
at 10~^ level performed with the Gravity Probe A spacecraft the deflection of light 
by solar gravity via very-long-baseline (radio) interferometry at 10~^ level [5]; the time- 
delay by gravitational potential using the Viking spacecrafts at Mars [6] and the Cassini 
mission at Saturn (the latter at 10~^ level ff]); the equivalence principle at 10"^'^ level and 
the geodetic precession at about six parts in 10"'^ with lunar laser ranging (LLR) of the 
Apollo and Lunokhod retroreflectors [8j; frame dragging with satellite laser ranging (SLR) 
of the LAGEOS satellites and with the Gravity Probe B (GPB) mission (the former at 10~^ 
level [9j); the latter is also expected to yield a very accurate measurement of the geodetic 
precession. 

The first three measurements above (perihelion precession, gravitational redshift and light 
deflection) are the three classical tests originally proposed by Einstein to verify his theory. 
The geodetic precession is a relativistic three-body effect that was predicted |10| by de Sitter 
in 1916 and observed in 1988 by LI. Shapiro et al. with an accuracy of 2% using LLR data 

m- 

In the effort of improving the experimental measurements and of possibly discovering new 
physics, several extensions and modifications of GR have been developed. One notable 
attempt is the Parametrized Post-Newtonian (PPN) formalism, whose verification has been 
the object of continuous and always improving experimental tests. No deviation from GR 
has been found so far. However, in presence of new physics beyond GR, it is natural to 
expect and try to measure modifications of the solar system observables described above, 
which historically marked the transition from newtonian physics to relativistic gravity. 
In this paper we treat the modification of GR to include spacetime torsion and motion 
along autoparallel orbits. Then we use the measured Moon geodetic precession and Mer- 
cury's perihelion advance to test our predictions. In a companion paper we show how 
the constraints on torsion provided by the measurement of the geodetic precession can be 
used also to constrain spacetime torsion with the frame dragging experiments on LAGEOS 
satellites. 

In the near future, improvements of the limits reported in these two papers may be ob- 
tained from the analysis of Mercury radar ranging (MRR) data taken since 1990 [1] , from 
the relentless analysis of more LLR data from more ground stations, from the mm-level 
range precision provided by the new APOLLO station at Apache Point, USA [13J (oper- 
ational since 2007) and from the release of the geodetic precession measurement by the 
GPB Collaboration. In the mid-term we expect a substantial advance from LLR with the 
deployment of 2""^ generation laser retrorefiector payloads with robotic soft-landings on the 
Moon, like the missions of the International Lunar Network (ILN), or similar geophysical 
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network, and like JAXA's Selene-2. We also expect that in the long term, more stringent 
limits can be set with the approved BepiColombo Mercury orbiter, an ESA Cornerstone 
mission. 

2 Theoretical framework 

An interesting generalization of GR includes a non-vanishing torsion. A class of theories 
allowing the presence of torsion is based on the extension of Riemann spacetime to Riemann- 
Cartan spacetime. The latter has a richer geometric structure, since it is endowed with a 
metric g^i/ and a connection F^^j, which is not the Levi-Civita connection. A compatibility 
condition between (7^,^ and T^^^, is required, namely the covariant derivative of the metric 
tensor must vanish identically. Under this assumption the resulting connection turns out 
to be non-symmetric, and such a lack of symmetry gives origin to a non- vanishing torsion 
tensor. We refer to |14j . |15| for the details. 

In most torsion theories of gravity, the source of torsion is the intrinsic spin of matter [TJ] , 
|16| . |17j . A recent review on searches for the role of spin and polarization in gravity can 
be found in [T8]. Since the spins of elementary particles in macroscopic matter are usually 
oriented in a random way, such theories predict a negligible amount of torsion generated by 
massive bodies. As a consequence spacetime torsion would be observationally negligible in 
the solar system. 

However, in |19j Mao, Tegmark, Guth and Cabi (MTGC) argue that, if there are theories 
giving rise to detectable torsion in the solar system, they should be tested experimentally. 
For this reason, in [19j a theory-independent framework based on symmetry arguments 
is developed, and it is determined by a set of parameters describing the torsion and the 
metric. Here, by theory-independent framework, we mean the following: the metric and the 
connection are parametrized, around a massive body, with the help of symmetry arguments, 
without reference to a torsion model based on a specific Lagrangian (or even on specific 
field equations). 

This framework can be used to constrain the above mentioned parameters from solar system 
experiments. In particular, MTGC suggest that GPB is an appropriate experiment for this 
task. In [19] the authors compute the precession of gyroscopes and find the constraints 
that GPB is able to place on the torsion parameters. In [20] Hehl and Obukhov argue that 
measuring torsion requires intrinsic spin, and criticize the approach of MTGC, since GPB 
gyroscopes do not carry uncompensated elementary particle spin. 

MTGC address also the question whether there exists a specific gravitational Lagrangian 
which yields a torsion signal detectable by the GPB experiment. As an example they quote 
the theory by Hayashi and Shirafuji (HS) in [21] where a massive body generates a torsion 
field. In such a theory gravitational forces are due entirely to spacetime torsion and not to 
curvature. The same property is shared by teleparallel theories [22], [23], [23], [25], [26] . 
Then MTGC propose what they call the Einstein-Hayashi-Shirafuji (EHS) Lagrangian, 
which is a linear interpolation of GR and HS Lagrangians. The main feature of the EHS 
theory is that it admits both curvature and torsion. 

The EHS model has been criticized by various authors. Flanagan and Rosenthal show in [27j 
that the linearized EHS theory becomes consistent only if the parameters in the Lagrangian 
satisfy suitable relations that, in turn, make the predictions coinciding with those of GR. In 
the paper [28] , Puetzfeld and Obukhov derive the equations of motion in the framework of 
metric-affine gravity theories, which includes the HS theory, and show that only test bodies 
with microstructure (such as spin) can couple to torsion. The conclusion is that the EHS 
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theory does not yield a torsion signal detectable for GPB. For this reason, in [19] the EHS 
Lagrangian is proposed just as a pedagogical toy model. In the present paper we will not 
treat the EHS model. 

As also remarked by Flanagan and Rosenthal in |27] . the failure of constructing the spe- 
cific EHS Lagrangian does not rule out the possibility that there may exist other torsion 
theories which could be usefully constrained by solar system experiments. Such torsion 
models should fit in the above mentioned theory-independent framework, similarly to a 
parametrized post-Newtonian framework including torsion. We remark that the parametrized 
formalism of MTGC does not take into account the intrinsic spin of matter as a possible 
source of torsion, and in this sense it cannot be a general torsion framework. However, it 
is adequate for the description of torsion around macroscopic massive bodies in the solar 
system, such as planets, being the intrinsic spin negligible when averaged over such bodies. 
For this reason we think it is worthwhile to continue the investigation of observable effects in 
the solar system of nonstandard torsion models within the MTGC parametrized formalism, 
under suitable working assumptions. In particular, our aim is to extend the GPB gyroscopes 
computations made in |19j to the case of motion of planets and satellites. 
In the present paper we compute, as an effect of spacetime torsion, the corrections both to 
the precession of the pericenter of a body orbiting around a central mass, and to the orbital 
geodetic precession. We describe the torsion by means of three parameters ti,t2,ts. Our 
computations show that a complete account of the precessions requires a parametrization 
of torsion up to an approximation order higher than the one considered in |19j . 
We consider the motion of a test body in a spherically symmetric field, and the motion 
of a satellite (either the Moon or LAGEOS) in the gravitational field of the Sun and the 
Earth. Since we use a parametrized framework without specifying the coupling of torsion to 
matter, we cannot derive the equations of motion of test bodies from the gravitational field 
equations. Therefore, in order to compute effects of torsion on the orbits of planets and 
satellites, we will work out the implications of the assumption that the trajectory of a test 
body is either an autoparallel curve or a geodesic. Such trajectories do not need to coincide 
when torsion is present. The computations will be carried out under the assumption of 
weak field and slow motion of the test body. 

We will assume that the motion of the satellite is obtained by superimposing the fields 
of the Sun and the Earth, both computed as if these bodies were at rest. Observe that 
these assumptions are satisfied to a sufficient order of approximation in classical General 
Relativity. 

As in the original paper of de Sitter [10], we characterize the motion using the orbital 
elements of the osculating ellipse. In terms of these orbital elements, the equations of motion 
then reduce to the Lagrangian planetary equations. We calculate the secular variations of 
the longitude of the node and of the longitude cu of the pericenter of the satellite. The 
computed secular variations show how the corrections to the orbital de Sitter effect depend 
on the torsion parameters , t2 > ^3 • In addition we calculate the secular variation of the 
longitude of the pericenter of a body orbiting around a central mass, and also in this case 
we find the corresponding dependence on ti, t2, t^. The data from the LLR and MRR 
measurements are then used to constrain the relevant linear combinations of the torsion 
parameters. 

Eventually, we consider the geodesic trajectories, and we find that torsion parameters cannot 
be constrained by solar system experiments. 

The paper is organized as follows. In Section [3] we briefly recall the notion of spacetime 
with torsion. In Section [4] we recall from [T9J how to parametrize the metric and torsion 
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tensors under the assumption of spherical symmetry, and we extend the parametrization 
up to a higher order of approximation. In Section [5] the connection up to the required 
order is given. In Section [6] we analyze the equations of autoparallel trajectories and we 
derive the related system of ordinary differential equations to second order. The expression 
of the system clearly reveals the perturbation due to torsion with respect to the de Sitter 
equations. In Section [7] we calculate the correction due to torsion to the precession of 
pericenter. In Section [8] we calculate the correction due to torsion of the third Kepler's law. 
In Section [9] we investigate the motion of a satellite in the gravitational field of the Sun 
and the Earth and we compute what we can call the perturbative forces due to torsion. 
In Section [TO] we derive the time evolution of the orbital elements of the satellite, using 
the classical perturbation theory of Celestial Mechanics, particularly the Gauss form of 
the Lagrange planetary equations. In Section 11 we calculate the secular variations of 
the orbital elements of the satellite. In Section 12 we give multiplicative torsion biases 
relative to the GR predictions. In Section 13 we report the constraints on the parameters 
of our torsion model from LLR and MRR, which is one of the main goals of this paper. In 
Section 14 we analyze shortly the geodesic trajectories. In Section 15 we summarize the 
future prospects of the LLR and MRR measurements and we discuss the implications of 
proposed and approved space missions for the search reported in this paper. Eventually, 
in the Appendix we confirm using a different formalism the computation leading to (7.11), 
and we show that, in the autoparallel scheme, torsion produces an effect on the precession 
of pericenter which was not taken into account in |60| . 



3 Spacetime with torsion 

We briefly recall the basic notions of Riemann-Cartan spacetimes [2], |15j . A spacetime 
equipped with a Lorentzian metric g'^i/ and a connection F^^^ compatible with the metric 
is called a Riemann-Cartan spacetime. Compatibility means that V ^g^jx = 0, where V 
denotes the covariant derivative. We recall in particular that for any vector field w'*' 

The connection is determined uniquely by g^^ and by the torsion tensor 



as follows: 



A 



rV = |^,|-^./' (3-1) 



where {•} is the Levi-Civita connection, and 



^y ~ ^uu^ (3-2) 



is the contortion tensor. In the particular case when L'^'^^, is symmetric with respect to /x, u 
the torsion tensor vanishes. In the present paper we will consider the case of nonsymmetric 
connections L^^^^. 



The Riemann tensor of the connection (3.1) is given by 



pu/j, ~ fj.p f^/^r + r j^Q,r r ^^.r (3.3) 

The particular case of vanishing torsion tensor corresponds to Riemann spacetime of GR, 
while the particular case of vanishing Riemann tensor corresponds to the Weitzenbock space 
time im. 
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4 Parametrizations of metric and torsion in spherical sym- 
metry 

Throughout this paper we use the natural gravitational units where c = 1 and G = 1. In 
the following we consider a spherically symmetric body of mass m. Introducing spherical 
coordinates {r,9,(p), we parametrize the metric and torsion tensors in a region of space 
(out of the body) where the dimensionless quantity Cm = m/r « 1 (i.e., r is large in 
comparison with the Schwarzschild radius of the body). As it will be shown in the sequel, 
such an approximation is accurate enough for the purpose of our computation^ 
We recall that Parametrized Post-Newtonian (PPN) calculations [29] show that a complete 
account of the pericenter precession must involve second order parameters in 6772 (for instance 
the PPN parameter /3). Therefore, assuming spherical symmetry, we parametrize the metric 
tensor ^^jy to second order. Under the assumption of spherical symmetry the line element 
has the following general expression in spherical coordinates: 

ds'^ = -h{r)dt^ + f{r)dr'^ + a{r)r'^[de^ + sin^ ed(l)\ (4.1) 

We can choose a{r) = 1, and to second order in m/r we have 

where T-L^J-^X^C are dimensionless parameters. The metric is then expressed in nonisotropic 
spherical coordinates. In the computations of trajectories that we will make in the following, 
only the function h{r) is required to the second order in Cm, while for /(r) the first order 
approximation is sufficient. 

We follow the notation of |19] for the parametrization of the metric tensor: the parametriza- 



tion (4.2) reduces to the first order in 6777, to tli6 one usGd. in ^19^ . In th.6 case of 3- PPN 



metric we have {30 ji Section 3.4.1]: 

^ = -2, J- = 27, X = 2(/3-7). (4.3) 



In the present paper all the others PPN parameters [29] are assumed to be negligible. 
When spacetime torsion is present, our calculations show that a complete account of the 



precessions must involve a parametrization of the torsion tensor Sau ^ up to second order in 



We now follow the spherical symmetry arguments used by [19]: to do this, it is convenient 
to parametrize the nonvanishing components of the torsion tensor in isotropic rectangular 
coordinates (t, ). We have 



s^k' = Qi: 



{1,2,3}, (4.4) 



where r' = (x^)^ + (x^)^ + (x^)^, V{r'), Q{r') are arbitrary dimensionless functions, and 
Sij is the Kronecker's symbol. If we consider such functions as depending on the dimension- 
less small quantity m/r' , for the purposes of our computations (see Section [5]) it is sufficient 

^For example, considering the field of the Sun of mass m and the Earth and Mercury as test bodies at a 
distance r of the order of their orbit semi-major axes, one gets, respectively, ~ 2 x 10~* and £,„ ~ 5 x 10~*. 
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to Taylor expand them up to the second order, hence we write 



m 



m 



m 



(4.5) 



where ti,t2,t3,i4 are dimensionless constants. In the particular case of first order approxi- 
mation, the above formulas yield the parametrization used in |19j . 

In order to transform (4.4) to nonisotropic spherical coordinates in which the metric (4.1) 
is expressed, it is convenient to first transform to isotropic spherical coordinates {r',9,(f)). 
We have for the non vanishing components of the torsion tensor 



— ti 

S„/n^ =S, 



m 



+ t3 
t2 



m 



m 



2(r')^ 



+ U- 



m 



(4.6) 



We now further transform (4.6) to nonisotropic spherical coordinates. To the required 
second order of accuracy, the transformation takes the form {t,r',9,(p) — )• {t,r,6,(j)) with 
r ~ r'(l + ^^)- The resulting expression of the components (4.6) of the torsion tensor in 
such coordinates is: 



J m m 



re 



m m 

2 



(4.7) 



The constants and t4 are related to ti, t2, and as follows: 



ts =t3 — —tl, 



~ T 

ti =t4 — — 12- 



(4.8) 



Therefore from (3.1) and (3.2) it follows that F^^^ becomes an explicit function of ti, t2, t^, 
ti, and the remaining four parameters involved. 



*2, ^3, ti, T-L, T ,T, C, r, 9, 



Since the metric and the torsion are constructed so that the compatibility condition 'V ^g\p 
is satisfied, then the metric parameters are independent of the torsion parameters. 
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5 The connection up to second order 



Using (3.1) and (3.2), the nonvanishing components of the connection in spherical symmetry, 
approximated to second order in €m — Tfijv^ read as follows: 



tr 



rt 



e 

i 4,, 
p9 

-r 

r' 



Br 



Tim 

'2^ 



n 

+ 



+ 



Tm 



2r' 



+ 



m 
2 



m ^ 
2 



2 



X + 2t3 



X 



m 



C 



m 



-r + {F + t2)m - {F^ + t2F + C- 2U) 



m 



-r sin^ e + {F + t2)m sin^ 9 - (F"^ + t2F + C 



\ 2 

2t4j — sm I 



r<f> 



t2- 



m 



2U 



m 



-sin^cos^, 
A cos 9 



smi 



In the computations of trajectories that we will make in the sequel only the components 
r*"^^ r*^^, r*^j are required to the second order in em, while for the remaining components 
the first order approximation is sufficient. The second order of approximation is used in 
Section [6] for F^'^j, and in the Appendix for r*^^, r*^j. It follows that the parameters C and 
(differently from X and t^) will not appear in the next sections, and consequently they 
will not be considered in the sequel of the paper. 



6 Equations of autoparallel trajectories 

The precise form of the equations of motion of bodies in the gravitational field depends 
on the way the matter couples to the metric and the torsion in the Lagrangian (or in the 
gravitational field equations). Here we consider the parametrized framework of Section 
|4] without specifying a coupling of torsion to matter, hence without specifying the field 
equations. 

In a Riemann-Cartan spacetime there are two different classes of curves, autoparallel and 
geodesic curves, respectively, which reduce to the geodesies of Riemann spacetime when 
torsion is zero [H]. Autoparallels are curves along which the velocity vector is transported 
parallel to itself by the connection T^^^. Geodesies are curves which are extremals of the 
length functional, and along which the velocity vector is transported parallel to itself by 
the Levi-Civita connection. In GR the two types of trajectories coincide while, in general, 
they may differ in the presence of torsion. They are identical when the torsion is totally 
skew-symmetric [14J , a special condition which is not satisfied within our parametrization. 
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The equations of motion of bodies in the gravitational field follow from the field equations, 
as a consequence of the Bianchi identities. The method of Papapetrou [31] can be used to 
derive the equations of motion of a test body having internal structure, such as for instance 
a small extended object that may have either rotational angular momentum or net spin. In 
standard torsion theories the trajectories of test bodies with internal structure, in general, 
are neither autoparallels nor geodesies [H], [15], [32], while structureless test bodies, such 
as spinless test particles, follow geodesic trajectories. 

In our computations of orbits of either a planet or a satellite (considered as a test body), 
we will neglect its internal structure. In a theory-independent framework we cannot derive 
the equations of motion from the gravitational field equations. Therefore we need some 
working assumptions on the trajectories of structureless test bodies: we will investigate the 
consequences of the assumption that the trajectories are either autoparallels or geodesies. 
Assuming the trajectory to be a geodesic is natural and consistent with standard torsion 



theories. However, we will see in Section 14 that the geodesies are the same as in the PPN 
framework. Hence new predictions related to torsion may arise only when considering the 
autoparallel trajectories, which will turn out to explicitely depend on torsion parameters. 
In the following we give some motivations which make worthwhile the investigation of 
autoparallel trajectories. 

In the paper [33], Kleinert and Pelster argue that the closure failure of parallelograms in 
the presence of torsion adds an additional term to the geodesies which causes structureless 
test bodies to follow autoparallel trajectories. Kleinert et al. also argue in [33], [3l] that 
autoparallel trajectories are consistent with the principle of inertia. Hehl and Obukhov in 
|20| criticized the approach of Kleinert et al., since the equations of autoparallel trajectories 
have not been derived from the energy-momentum conservation laws. Kleinert investigates 
this issue in |35j where the autoparallel trajectories are derived from the gravitational 
field equations via the Bianchi identities, in the case when torsion is derived from a scalar 
potential (see [15] for a discussion of such a kind of torsion) . 

In the papers [36], [37], using a reformulation of the theory of Brans-Dicke in terms of a 
connection with torsion [3E], Dereli, Tucker et al. suggest that the autoparallel trajectory 
of a spinless test particle is a possibility that has to be taken into account. In |36j the 
results of the investigation of autoparallel trajectories are applied to the computation of the 
orbit of Mercury. In [39], [40j the equations of autoparallel trajectories are derived from the 
gravitational field equations and Bianchi identities, in the special case of matter modeled as 
a pressureless fluid, and torsion expressed solely in terms of the gradient of the Brans-Dicke 
scalar field. 

The above mentioned results show that there is an interest in the autoparallels in spacetime 
with torsion, which make worthwhile their investigation in the present paper. The system 
of equations of autoparallel trajectories of a test body reads as 

where r is the proper time |41j. Notice that only the symmetric part ^(r'*' + T'*' ) of the 



measured from (6.1) 



connection enters in (6.1); in addition, the totally antisymmetric part of Sx^^y cannot be 



The trajectory of a test body has to be a time-like curve. Since the connection is compati- 
ble with the metric the quantity g^y^^^^ is conserved by parallel transport. The tangent 
vector to the trajectory undergoes parallel transport by the connection along the au- 
toparallel. Therefore, an autoparallel that is time-like at one point has this same orientation 
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everywhere, so that the trajectory is strictly contained in the hght cone determined by g^^, 
in a neighbourhood of every of its points. Hence the compatibihty of the connection with 
the metric ensures that autoparallels fulfil a necessary requirement for causality. 



The equations (6.1) can be rewritten as 



df^ y f"" ^ f"" dt J dt dt 

for a e {1, 2, 3}. In our units (^)^ and are of the order of Cm- 



(6.2) 



We use for x° spherical coordinates {r,9,(f)). Substituting in (6.2) the expression of F 



given in Section [s] one gets, to the order of accuracy. 



m 



, m 



,m 



.m 



+ + + C^r^ + (r + Vm) 6*^ + (r + Vm) sin^ 6 



2 ^771 



2 „m 



rO + sin ^ cos cj)^ , 
i-cf)- 2 cot 6 9 4), 



(6.3) 



where 



A 



B 



+ Z-ti{U-F)-2^, 



n + 



T 



V 



£ = ti + t2-n. 



In order to take into account relativistic corrections, the right hand sides of (6.3) must be 



at least of second order in em', this is guaranteed if F^^^ is developed to second order in em, 
while it is enough to develop the remaining components of the connection to the first order. 



System (6.3) to lowest order becomes 



dv 
dt 



ti 



n 



m 



(6.4) 



where ir is the unit vector in the radial direction. Imposing the Newtonian limit it follows 
(see also [IHl formula (23)]) 

h-^ = l, (6.5) 



hence ^ = 0. 
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We now transform (6.3) in rectangular coordinates x = r sin 9 cos (j), y = rsin^sini?^, z 
r cos 9. Writing x^,x''^, for x, y, z we get 



X 



-m^ + Dm —r+cm — ^+I'm^q? + -J^m^r" 



a = 1, 2, 3, 



(6.6) 



Note that in case of no torsion (i.e., ti = t2 = ^3 = 0) and when T = 2 and X = 
system (6.6) reduces to the equations of motion of General Relativity in the weak field 
approximation. 



7 Precession of pericenter 

From the second equation in (|6.3|) it follows that if 9 and 9 vanish at one time then 9 is 



identically zero. Therefore, assuming plane motion, the system (6.6) can be written in the 
form 



m 



m 



X -\- Fxi 



y 



J/' 



where {Fx,Fy) is the perturbation with respect to the Newton force, 

F, 



Jut T 2 -I— 2 

15mf —r ^ cm^ + um^^ H — J-m^r , 

fj^^ rjr^^J 2 y^O 



Bm 



(7.1) 



and now = + y^. 

The vector (F^,, F^) can be decomposed in the standard way along two mutually orthogonal 
axes as 



S — -Fx + -Fy, 
r r 



T 



dix/r) ^ diy/r) ^ 



(7.2) 



du ^ ' du 

Here S is the component along the instantaneous radius vector, T is the component per- 
pendicular to the instantaneous radius vector in the direction of motion, where u is the 
argument of latitude. Then, substituting (7.1) into (7.2) gives 

,2 ^2 



S = + Cm-^ + Vmii 



(7.3) 



ru 



Let us now recall 



that, using the method of variation of constants, 

a(l -e2) 

r = , 

1 + e cos V 



(7.4) 
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where a is the semimajor axis of the satelhte orbit, e is the eccentricity, v is the true 
anomaly, and 



r e smv 



2- 2/1 2\l/2 

r V = na [1 — e ) ' , 



(7.5) 



a(l-e2) ' 

n = 2tt/U, U the period of revolution. Recall that, in the Newtonian approximation, 
= mjc? from the Kepler's third law. Following the standard astronomical notation, 
we let bj be the longitude of the pericenter, and u = v ^Zj. We also recall the following 
planetary equation of Lagrange in the Gauss form [43, Ch. 6, Sec. 6]: 



diij 
~dt 



(l-e2)i/2 



nae 



-S cos V + T { I + 



a(l-e2) 



smv 



(7.6) 



Notice that, since 5 and T are of the order e^, we have that ^ is of the order Cm^- We are 
therefore allowed to make the approximation 



(7.7) 



Inserting (7.3) into (7.6), making also use of (7.4), (7.5), (7.7) and the Kepler's third law, 
we have 



duj 
~dt 



Bm? Cem n 

-(1 + e cost") cos ft; tz ftt sm vcosvv 



a(l-e2) 



(7.8) 



v?a'^e{l — e^) 

T)m 2 "^"^ 2 / 

^(1 + ecosu) cosTJ-uH — ^ sin v{2 + ecosv)v. (7.9) 

ae(l — e^) o(l — e^) 

According to perturbation theory, we now regard the orbital elements on the right hand 



side of (7.8) as approximately constants. Therefore, integrating with respect to t, we obtain 
Bm 



6u} 



( e e 

ae{l - e^) V 2 4 



Cem 



3o(l-e2; 



sm V 



1 



— ^ sin f + H — sin(2f ) + e sin v sm" v 

ae(l — e^j \ 2 3 



a 1 



V sin(2f ) H — sin v . 



(7.10) 



7.1 Correction to the precession of pericenter due to torsion 

Secular terms appear in buj. Using the expressions for T) and such secular contributions 
are: 



If ^1=^2 = ^3 = 0, using (4.3) then we find 

(55)sec = (2 + 27 - /3) 



a(l 



(7.11) 



(7.12) 



which gives the precession of pericenter in terms of PPN parameters, as it can be found in 
[291 Chapter 7, formula (7.54)]. Moreover, when = -2, J" = 2 and X = (i.e., ^ = 7 = 1) 
we find the usual expression of ((Jw)^^ given by General Relativity. In the case of Mercury, 
such a precession amounts to 42.98 arcsec/century. 
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Our formula for 5u}^pr differs from the formula 



(5w 



,GR 
'sec ' 



(7.13) 



found by Mao et al. in [19, formula (CIO)]. Formula (7.13) does not reproduce the PPN 
result (7.12) when the torsion parameters vanish, though it reproduces the GR result in 

-2 and T = 2. In the Appendix we compute the precession of 



the particular case 'K 
pericenter following the method used in [T9], obtaining again expression (7.11). 
The precession of the pericenter in a reformulation of the Brans-Dicke theory in terms of a 
connection with torsion has also been computed in [36j by using autoparallel trajectories. 



8 Correction to Kepler's third law 

In this section we compute the relativistic correction of Kepler's third law for Earth motion, 
in the presence of torsion. This result will be used in the sequel, for the computation of 
the satellite geodetic precession. We note that such a correction was not necessary in the 
previous computation of the precession of pericenter at the required order of accuracy. 
We introduce the following coordinates. We take a system of rectangular coordinates cen- 
tered at the Sun. The triplet (^, r/, Q denotes the Earth's coordinates in this system, and 
we assume that the ecliptic plane coincides with the plane = 0; we assume that the 
eccentricity is zero, so that the Earth's orbit is given by 



^ = pcosL, r] = psmL, C = 0) 
where p and L denote the radius of the orbit and the longitude of the Earth, respectively. 



Therefore, system (6.3), using 9 = tt/2 and p = 0, yields 

- -^+B^ + {p + Vm)L^ 



0, 



i.l) 



where m, here and in the sequel, denotes the mass of the Sun, supposed spherically sym- 
metric. It follows that L = uq, where uq is constant, and that 



m 
7^3 



l-{T + t2) 



U. 

2 



l-{J' + t2) 



12) 



where we have used the Newtonian limit condition (6.5). Since the semi-major axis of 



the Earth's orbit is large with respect to the Schwarzschild radius of the Sun, we have 
m/p ~ 2 X 10~^ << 1. Since in our units c = 1, we also have that Vq « 1. It follows that. 



up to second order, equation (8.2) becomes 



m 
"~3 



2 



1 



){F + t2)-B 



m 



2 



p 



-I 



l-{F + t2-B) 



m 



P 



13) 



This approximation will be used for the computation of the satellite geodetic precession in 
the next section. 
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9 Motion of a satellite in the field of the Sun and the Earth 



In this section we investigate the motion of a satelhte (either the Moon or LAGEOS) in the 
gravitational field of the Sun and the Earth in presence of torsion. The coordinates ry, 
have been defined in Section [sj The triplet {X, Y, Z) denotes the satellite's coordinates, and 
we set 

a2 = + + ^2. 

The satellite's coordinates with respect to the Sun will be written as 

X = C + x, Y = ri + y, Z = ( + z, 

where {x, y, z) are the coordinates of the satellite with respect to the Earth. We use the 
standard coordinates transformation 1421. 1431 used in Celestial Mechanics 



X = r (cos u cos Q — sin u sin cos i) , 

y = r(costisinil + sinucosr^cosi), (9-1) 
^ z = r sinwsinz, 

where r is the distance between the Earth and the satellite, u is the argument of the latitude, 
fl is the longitude of the node, and i is the orbital inclination. 

We suppose that the Earth is spherically symmetric. The semi-major axes of the Moon 
and LAGEOS orbits around the Earth are small in comparison with p, so that in our 
computations we will neglect the powers of r/p greater than on^ Hence, also m/A << 1. 
We will assume that the motion of the satellite is obtained by superimposing the fields (i.e., 
the connections as described in Section [5| of the Sun and the Earth, both computed as if 
these bodies were at rest. More precisely we assume 

9fiu = ig,iu)o + {gfiu)l S^j,^ = {S^^^)o + iS^j,^)i, (9.2) 

where 

- {gfiu)o and (5'^j^'^)o are the metric and the torsion tensors, as given in Section|4| taking 
into account the Sun only, supposed at rest; 

- {9fiu)i and {S^^^)^ are the metric and the torsion tensors taking into account the 
Earth only; these tensors are given at each time by the expressions in Section |4j 
computed as if the Earth were at rest (at that time). 

Observe that these assumptions are satisfied to a sufficient order of approximation in clas- 
sical General Relativity [lOj, [44J if all the terms that give rise to periodic perturbations are 
neglected (we are interested only in secular effects). We also note that a rigorous justifica- 



tion of the validity of (9.2) would probably require an extension of the parametrized torsion 
model to the case of three interacting bodies, which is beyond the scope of the present 
paper. 

As a consequence we have 



d^_f(Px^y / d2X" \ _(d^ 



where we write ^ for X, y, z. Similarly we write X^,X'^,X^ for X, Y, Z, and 4^, 4'^ 

for ^, C- Moreover 



^For the Moon and LAGEOS we have r/p ~ 2.6 x lO""^ and r/p ~ 8.5 x lO""^ respectively. 
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^® ^^^^ hand side of (6.2) with the coefficients of the connection computed 



using {g^^)1 and {S ■ 



( ) hand side of (6.2) with the coefficients computed using {g^u)o and 



- ^^^^ is the left hand side of (6.2) with the coefficients computed using {g^u)o and 



The contribution of the term ^^^^ gives a secular precession of the perigee which has 

been computed in Section [7} The other terms represent the perturbing effect of the Sun. 
Since all the perturbations here considered are small enough to allow us to superimpose 
them linearly, in what follows we compute the perturbing effect of the Sun only. With these 



assumptions, using the Newtonian limit (6.5) the right hand members of (6.6) give 



x° + = kAA'^ + ksB" + kcC" + koD'^ for a = 1, 2, 3, 



(9.3) 



where 



A" 



m 



m 







A3 " 








" 






and for notational convenience we set 

kA = B, ks = £, kc = V, ko 



37" 



(9.4) 



The left hand side of (9.3) contains the ordinary Newtonian perturbing function N"^, which 
also requires a correction , according to the computations in Section [sj This will be 
made clear in Section [lOl 

The components R°' and j4°, approximated to the first order with respeclj^to r/p and r/p, 
read as follows (where we write Rx , Ry and Rz for R^ , R^ and R^ respectively, and similarly 



^From formula ( |10.10[ ) below expressing r, it follows |r| < -^^j so that r is at least small as w, hence \r/p\ 
is smaller than r/p. 
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for A,B etc.): 



Rx = ^{x- Sxcos^ L) + Pr^, 



Ry = 

Rz = 

■^x — 

Ay = 

A., = 



m 

7 

m 

~7 

2 

m 



{y-Sysin^L) + Pr^, 

(a;-4xcos^ L) + Pa^, 
[y-Aysin^L] + P^,, 
-z + Pa.. 



(9.5) 



(9.6) 



Here the terms P^a , Pj^a denote a finite sum of addenda with the fohowing property; each 
addendum is of the form / (sin L)"! (cos L)"^^ where ni + n2 is odd and / is a factor inde- 
pendent of L. Such terms will give periodic contributions to the perturbations of the orbital 
elements: therefore they will be neglected in the computation of secular perturbations. 
The components B", C" and D", approximated to the first order with respect to r/p, r/p, 
and taking into account also the terms in rf/p^, read as follows: 



Bx = 



mL 
P 

mL 



{Lx-y) s\T? L + Pb^, 
{Ly + x)cos^ L + Psy, 



(9.7) 



Cx 



Cz 



I Bz = Pb., 
mL 



P 

m,L 
P 



Lx + {-3Lx + 2y) cos^ + Pc, , 

Ly + i-3Ly - 2x) sin^ l) +Pn, 



(9.8) 



-z + Pc., 






' Dx 


= Pd. 


< 


Dy 


= PDy 






= Pd. 



(9.9) 



Also here the terms Pb" , Pc" and Pd" have the same structure of Pro and Pa" , thus giving 
periodic perturbations of the orbital elements. 
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10 Computation of orbital elements via perturbation theory 



In this section we introduce the tools from Celestial Mechanics needed to compute the 
secular perturbations of the orbital elements. In the following the orbital elements, the true 
anomaly v and the argument of latitude u will be referred to the satellite's orbit around 
Earth. 



Using the Newtonian limit condition (6.5) and the correction to Kepler's third law (8.3) we 
have 

i?" = AT" _ fe^TV" for a = 1,2,3, 



where 



' A^^ = u^x (1 - 3cos2 L) + Pn,, 



Ny = iy^y {1-3 cos^L) +Pn^, 



2 



(10.1) 



(10.2) 



P 



y(l-3sin2L)+P, 



(10.3) 



-z + P, 



and 



kN = + t2 - kA- 



(10.4) 



Again the terms P/v" and give periodic contributions to the perturbations of the orbital 
elements. 



In equation (10.1) is decomposed into the ordinary Newtonian perturbing function N'^ 



plus a relativistic correction —k^N", which depends also on torsion. 



Equations (9.3) can be rewritten as 

x'^ + Ar° = i?" for a = 1,2,3, (10.5) 
where is the perturbation with respect to the Newton force, and it is given by 



Recalling also (9.1), the perturbation {Fx, Fy, Fz) can be decomposed in the standard way 
along three mutually orthogonal axes as 



X ^ y ^ z ^ 
S = -Fx + -Fy + -F„ 



d{x/r) , d{y/r) ^ , d{z/r) ^ 



du 
sinu W = 



du 



du 



(10.6) 



d{x/r) d{y/r) d{z/r) 

Ol Ol Ol 



Here S is the component along the instantaneous radius vector, T is the component per- 
pendicular to the instantaneous radius vector in the direction of motion, and W is the 
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component normal to the osculating plane of the orbit (collinear with the angular momen- 
tum vector). 

Recalling from Section [s] that L = uq, L = uqI^ we replace S, and T with the averages 



27r 



Wdt, 



27r 



Sdt, 



2-K 



Tdt, 



respectively. Taking these averages has the following consequences: (i) it eliminates the 
dependence of W on the trigonometric functions of L, hence the periodic terms in (9.5)- 
(9.9) disappear; (ii) the remaining terms, contributing to the secular effects, are multiplied 
by a factor depending on vq. 

In order to compute the components S, T and W of the perturbation, using the method of 
variation of constants HH, 1431 we write 



X = r(cos u cos — sinnsinri cosi) + ru{— sin u cos — cos u sin cos i), 
y = r(cos usinri + sin it cos Q cosi) + ru{— sin n sin + cos u cos cos i), 
z = f sin u sin i + r cos u sin i, 

in order to eliminate i". 

We recall that u = v + u} — (where w is the longitude of the pericenter) , which allows to 
make again the approximation u = v, and we make use of the area law xy — yx = r'^vcosi 
in order to simplify the computations. 
We decompose 

W = Wa + Wb + Wc + Wf^, 

with obvious meaning of the notation. 
We have 

' 2 

TTl 



WB = kB 

Wc = kc 



P 



J mvl 3 

Wf, = k^ — - 



{—UQZCosi + VZ — r cosusmi) 



cos i — vz + f cos u sin i 



zcosi. 



Rearranging terms it follows 



W 



p 



ks , 3% , 'ikN\ 

Ik A z cos I 

2 2 2 



+ 



mfou ( k 
P 



kc]z + 



k 



B 



+ kc r cosusini. 



2 / p 

When ^ = -2, J" = 2 and T = (i.e., /? = 7 = 1) and ti = t2 = ^3 = 0, we find 

3mz^o 



Wgr = zv 

P 



-r costismz. 



P 



(10.7) 



When the satellite is the Moon, we have rsini = 0{ei) which is negligible. In this case 
W^GR corresponds to the formula found by de Sitter in [lOl (95)]. 

Now we compute the Gaussian component S of the perturbation. Analogously we decom- 
pose 

S = Sa + Sb + Sc + 'S'jy- 
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We have 



SA = -kA'^{r'-2z% 
rp 



S 



B = KB 



Sc = k 



2rp 



rp 



V COS I) — UqZ 



c- 



^0 



+ V COS ^ + -I'oz 



(10. 



Rearranging terms it follows 



rp 



ks kc kjsf \ 



ks 



+ kc ] vr cos i 



When n 



I ks 3 3 \ 2 

+ I 2kA - + -jkc + -kN ) i^qz 



-2, T = 2 and T = and ^1=^2 = ^3 = 0, we find 



dgr = —3 rcos^. 



When the satellite is the Moon if we approximate cosi ~ 1 then 5gr corresponds to the 
formula found in [10' (95)]. 
Similarly we decompose 

T = Ta + Tb + Tc + T^. 



We have 



m 

Ta = 2A;yi— j-zcosnsini, 

J-B = ks (—^qz cos Msmz + r cos i) 

2p 

mvQ ( 3 

J-C = kc -z^o-2cosusmz — rcos^ 

p \2 

3 mun 
N ~ 2 zcostismz. 



P 



Rearranging terms it follows 
rnvQ 



P 



kB 3 3 
2kA ^ + 2^"^ + 2 N I J^o-zcosusm? + 



kB_ 
2 



(10.9) 



kc rcosi 



When n 



-2, = 2 and 1 = and ^1 = ^2 = ^3 = 0, we find 

J-GR = 3 rcosi. 



When the satellite is the Moon if we approximate cosi ~ 1 then Tqr corresponds to the 
formula found in |lOi (95)]. 

As for the computation of the precession of pericenter in Section [7j we use the formulae 

"'^'^ . ae(l — e'^)vsmv 



a(l - e\ 
1 + e cos V ' 



(1 + e cos f )^ 



(10.10) 
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We also recall the following planetary equations of Lagrange in the Gauss form [43^ Ch. 6, 
Sec. 6]: 



dt 



1 



na^(l — e^)^/^ sini 
du (1 - e2)V2 



Wr sinn, 



flO.lll 



dt 



nae 



S cos V + T { 1 + 



ail 



smv 



+ 2 sin^ 



1 di} 

2 ~dt' 



where n = 27r/U , U the period of revolution of the satellite around Earth. 

We make the computations up to the first order with respect to the eccentricity e, and we 

use the following formula for the true anomaly |431 formula (2.6.7)] 

v{t) = n{t - r) + 2esin [n{t - r)] + 0(e), (10.12) 

where r is the satellite time of perigee passage. 



11 Precession of orbital elements in the presence of torsion 

Using the expressions of S, T and W computed in the previous section and integrating 
(10.11) we find secular terms in the expressions of 5^1 and 6(jJ. According to perturbation 
theory, we regard the orbital elements as approximately constant in the computation of 
such integrals. Since u = v + oj — fl, we can make use of the approximation iic^i v. Moreover 

p2U/2 



we use the formula na (1 — e 



rH, and (10.12). 



Let us first consider the computations for the node il. Substituting the decomposition 

dt 



(10.7) of W into the expression of ^ in (10.11), one has to compute the following three 



types of integrals: 
^ , zr sin u 



na2(l-e2)V2 



dt, 



II 



zvr sm u 



na?{l - e2)i/2 



dt. 



III 



rr sm u cos u 



na2(l-e2)V2 



dt. 



The integrals I and II yield periodic terms plus the following secular contributions: 



Sln^ 



n 



t 



V 

2n 



lis 



Sim 



.t 
2' 



where t is time. The integral III yields only periodic terms. 

In conclusion, it turns out that the secular contributions to the variation of the node are: 



P 



kB , Sfcc , 3kN\ 

2kA 

2 2 2 



t 



V 

2n 



cos I ^ mvQ 
n p 



2 



kc 



(11.1) 



Since v = nt + periodic terms in v, inserting (9.4) and (10.4) into (11.1) we obtain 

n 




2 2 2 



1^0 

cos ^ 

2n 



n 



{U + F)-Z-tiF + ti + t2 + 2tz 



(11.2) 



t. 



Using the Newtonian limit (6.5) and setting 

n 



Ci = 1 



+ 2.F + 3t2, 



Co 



, n -^2 



T -T + t2 + 2t^, 



(11.3) 
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we obtain 



When torsion is zero, that is ti 
(PPN formalism) we obtain 

1 rnvQ 



1 rnvQ 

t2 = h = 



G2 — cos I I t. 

n 



0, if we let J- = 27, X = 2(/3 - 7) and Ti 
1 mz^o n\ ^0 



(1 + 27) t - ^— (1 - /3) — cosi t. 
2 p n 



(11.4) 
= -2 

(11.5) 



The first term on the right hand side of (11.5) determines the usual geodetic precession 
effect. The second term is consistent with the computations in [IS, formula (48A)], when 
the satellite is the Moon, so that we can approximate cosi ~ 1. 

Letting 7 = /3 = 1 we find the usual formula of geodetic precession in General Relativity 
found by de Sitter (see [TOl formula (97)]), 



GR 
sec 



2p 



-t. 



We recall that when the satellite is the Moon, this precession amounts to 1.92 arcsec/century. 
Now we consider the computations for the perigee. The contribution of the gaussian com- 
ponent S to the variation 60J of the perigee is given by the integral 



(l-e2)i/2 



nae 



Scosv dt. 



(11.6) 



Substituting the decomposition (10.8) of S into the above integral, one has to compute the 
following three types of integrals: 



r cos V dt. 



rvcosv dt, 



cost" dt. 



We evaluate such integrals for small values of the eccentricity e and, taking into account 
that e appears at the denominator of (11.6), we expand the integral up to the second order 
in e. This is accomplished by expanding r in (10.10) with respect to e, and inserting (10.12) 
in the resulting expression. Then each of the integrands of I*^, II"^, III'^ turns out to be a 
sum of products of simple trigonometric functions, from which secular terms and periodic 
terms can be separated. Such integrals then yield the following secular contributions: 



r 



-ae(l — e^)t, 



ir 

-L-I-C, 



1 



ae(l — e^)f , 



and 



III'- 



8 



ae(l — e^) sin^ i [3 sin^ (oj — ^l) + cos^ {uj — Jl)] t. 



(11.7) 



In order to evaluate the integral III"^, since z = r sin u sin i, we have made use of the relation 
u = v + ui — to express the trigonometric functions of u. 

The contribution of the gaussian component T to the variation 6uj of the perigee is given 
by the integral 



(l-e2)i/2 



nae 



T 1 + 



a(l 



sinv dt. 



Substituting the decomposition (10.9) of T into the above integral, one has to compute the 
following two types of integrals: 



I^ 



1 + 



a(l-e2) 



zcosiisinw dt, 



1 + 



a(l-e2) 



rsinv dt. 



21 



As in the previous case, we evaluate such integrals for small values of the eccentricity e. 
Moreover, we will use the expression in (10.10) for r and the relation u = v + uj — Vt. 
Then each of the integrands of I"^, 11^ turns out to be a sum of products of simple trigono- 
metric functions, and we can extract the following secular contributions: 

Isec =-oe(l — e^) sini [sin^ (w — 17) — cos^ (w — 17)] t. 



II 



=ae(l — e^)v. 



In conclusion, the secular contribution to the variation of the perigee w is: 

ks kc kN \ 



1 mvQ / ks \ 3 mvl 
'2 p \2 ^ ) 2 np 



-kA + 



+ 



r5sin2isin2(S3-f7)-(l-cosi)] ( 2k 
2 np ^ \ 



ks 3, 3, , 



(11. 



Inserting (9.4) and (10.4) into (11.8) we obtain 

3 mvn 



1 mvQ ( % 



2 p 



2 +-^+2+ 2 



t + 



4 np 



-n- T-i + h + t2 + 2t3 + tiU) t 



1 mvQ 
4 np 



[5 sin^ i sin2(w - 17) - (1 - cos i)] (^{U + F)-!- tiF + ti + ^2 + 2t^ t. 



Using (6.5) and (11.3) we get 



4 p I n 
When torsion is zero, in the PPN formalism we obtain 



{Ci + C72^ [4-cosi-5sin2isin2(5-17)]} t. (11.9) 



1 + 27 l-/3^t 
2 p I p n 



1 mvl{ 



2 np 



[5 sin^ i w^ioj - 17) - (1 - cos i)] (1 - /3)t. 



(11.10) 



The first term on the right hand side of (11.10) determines the usual geodetic precession 
effect. The second term is consistent with the computations in |45t formula (48A)], when 
the satellite is the Moon, so that we can approximate cos i ~ 1 (and the third term becomes 
negligible). 

Letting 7 = /3 = 1 we find Ci = 6, C2 = 0, hence the usual formula of geodetic precession 
in General Relativity found by de Sitter (see [lOl formula (97)]), 



(5w 



>GR 
'sec 



SniVQ 

! 

2p 



Both in (517)sec and in {5uj)sec there appears the term 



1 muQ 



1 - ^ + 2^ + 3t2 ) i, 



which is independent of n and thus of the details of the satellite's motion. This term can 
therefore be interpreted as the geodetic precession effect when torsion is present. 
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12 Torsion biases 

In this section, similarly to [19j, we define multiplicative torsion biases relative to the GR 
predictions. These torsion biases will be used to put constraints on torsion parameters from 
solar system experiments. 

For the case of precession of the satellite orbital elements we define 



bo.= 



GR' 



GR- 

sec 



(12.1) 



From ( |11.4[ ) and ( |11.9[ ) it follows 

1 



bn =^ (Ci- C2— cost 
b \ n 



6~ =i |c7i + ^2 — [4 - cos i - 5 sin^ i sin^ (S - Vl)] \ 
o y. n '- -'J 



(12.2) 



where the constants Ci and C2 are defined in (11.3). 

For the purpose of comparison of solar system experiments with the predictions from the 
present torsion theory, we will assume that all metric parameters take the same form as in 
PPN formalism, according to (4.3). We have 



Ci = 2 + 47 + 3*2, 

C2 = 2 - 2/3 + t2 + 2*3, 



(12.3) 



and from (|12.2|) we obtain 
bn 



1 



3(1 + 27)+ 2 



COS Z 1- fj + — + ts 



3n 



1 . s *2 1^0 r 



5 sin^ i sin^ {u - n)] |^1 - /3 + ^ + *3^ . 



(12.4) 

Note that the torsion correction to the geodetic precession in (12.4) (namely the term 
^(1 + 27) + ^ ) differs from the corresponding one found for gyroscopes in \W\ formula 
(47)] by a numerical factor of order of unity. 

When the satellite is the Moon, we have cosi = 1 + O(i^), sin^i = O(i^), and we may 
approximate 



1 / X *2 t'n 

-1 + 27) + -- — 
3^ 2 3n 



t2 , T^O 



6s = -(1 + 27) + ^ + ^ 1 - /3 + ^ + *3 



n 



1 - /3 + I + t3 ) , 



*2 



(12.5) 



Similarly, considering now du as the precession of pericenter computed in Section [7| we 
have 



Br, = 



3 V'^^ 4 ^ 2 2 



In the case of a PPN metric we find 

1 



= - (2 + 27 - /3 + 2*2 + *3) ■ 



(12.7) 
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13 Constraining torsion with the Moon and Mercury 



Here we compare the predicted torsion biases to experimental measurements in order to set 
limits on the torsion parameters. 

Recent limits on various components of the torsion tensor, obtained in a different tor- 
sion model based on the fact that background torsion may violate effective local Lorentz 
invariance, have been obtained in [l6]. See also [57], where constraints on possible new 
spin-coupled interactions using a torsion pendulum are described. 



13.1 Moon: geodetic precession 

The GR test of the geodetic precession, evaluated with LLR data and expressed as a relative 
deviation from the value expected in GR, is Kgp = —0.0019 it 0.0064 [8]. In our torsion 
theory this is to be compared with the first two terms on the right hand side of equation 



( |12.4D : 

I - 1| = I J(l + 27) + I - 1| = 1^(7 - 1) + || < 0.0064, (13.1) 



where, taking into account the last sentence at the end of Section 11, we define 



fe = ^(l + 27) + |. 

Using the Cassini measurement 7 — 1 = (2. lit 2.3) x 10~^, we can neglect this term compared 
to the experimental uncertainty on Kgp and get the following constraint on t2- 

2 I 6^ - 1| = Ital < 0.0128. (13.2) 

The meaning of the constraint on the torsion parameter t2 is the following. Using the value 
((^a;)^^ = 1.92 arcsec/century for the geodetic precession of the Moon's orbit in GR, the 
geodetic precession in the presence of torsion is 



1.92 arcsec/century. 



If the parameter t2 had a value larger than 0.0128, this would imply the precession of the 
Moon's perigee would be (neglecting the contribution of 7 — 1) 

{6ui)sec > (1 + 0.0064)1.92 arcsec/century, 

which would be inconsistent with LLR data. If t2 < —0.0128 we would have an analogous 
inconsistency. To give the reader a further feeling of the effect of a nonzero value of t2 
in terms of orbit displacement, we provide the following, extreme example. The 1.92 arc- 
sec/century precession amounts to a perigee displacement of about 3 meters per lunar orbit 
period (about 27 days). A value t2 = 1 would imply (neglecting the contribution of 7 — 1) 
a geodetic precession in the presence of torsion of about 4.5 meters/lunar orbit period. 
The relentless accumulation of LLR data, the mm range accuracy of the APOLLO station, 
the start or restart of LLR operation of additional ILR^ stations (like MLRO, the Matera 
Laser Ranging Observatory in the south of Italy) will provide continuous further improve- 
ments of the Kgp test and therefore, of this limit on t2. Future improvements are possible 
also with current data and stations, by further developing and refining the current orbital 
software packages. 



International Laser Ranging Service; see http://ilrs.gsfc.nasa.gov/ 



24 



13.2 Mercury: perihelion advance 



The perihelion advance of Mercury has been measured with planetary radar ranging by 
Shapiro et al. in 1989 [3]. They found it to be consistent with its GR value with a relative 
standard error of 10~^. In the PPN framework, this can be used to "infer that /3 = 1 to 
within a standard error of (t(/3) = 0.003" (quoted from [3j). 

According to our torsion model, in the case of the PPN metric with the torsion bias given 



by (12.7), we get for Mercury the one standard deviation limit: 



l^ij - 1| = ^|27 - /3 - 1 + 2^2 + tsl < 0.001, 



|27-/3-l + 2t2 + t3| < 0.003, 



(13.3) 



1(27 - 2) + (1 - /?) + 2*2 + h\ < 0.003. 



Using the Cassini measurement 7 — 1 = (2.1ib2.3) x 10 ^, we can neglect this term compared 
to the experimental uncertainty on the perihelion advance and get the constraint: 



|1 -/3 + 2t2 + t3| < 0.003. 



(13.4) 



The limits (13.2) and (13.4) on the values of t2 and (1 — (3) + t3 are represented graphically 



m 



Combining the LLR and MRR constraints one gets on the following limit: 

|1 -/3 + t3| < 0.0286. 



(13.5) 



If we assume the Nordtvedt effect |48) . and that the Nordtvedt parameter rjj^ = 4/3 — 7 — 3, 
then the measured value [8j is /3 = 1 + (1.2 it 1.1) x 10~^, which makes (3 — 1 negligible 
compared to the experimental uncertainty on the perihelion advance. The constraint then 
becomes: 

|2t2 + ^31 < 0.003. (13.6) 

In this latter case, combining the LLR and MRR constraints one gets on the following 
limit: 

|t3| < 0.0286. (13.7) 

The meaning of the constraints on the torsion parameter ^3 is the following. Using the 
value (5tD)2^ = 42.98 arcsec/century for the precession of Mercury's perihelion in GR, the 
precession of the perihelion in the presence of torsion is (neglecting the contribution of 
7-1): 

1 



(5(5) 



1 + - (1 - /3 + 2t2 + t3) 



42.98 arcsec/century. 



If the linear combination (1 — /3 + 2*2 + ^3) had a value larger than 0.003, this would imply 
the precession of Mercury's perihelion would be 



i5u))sec > (1 + 0.001)42.98 arcsec/century. 



(13. 



which would be inconsistent with MRR data. If the parameter t2 takes the lowest value 
consistent with LLR data, i.e., if t2 = —0.0128, then a value of the parameter larger 
than (0.0286 + /3 — 1) would imply the inconsistency (13.8) with MRR data. If /3 — 1 is 
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-0.03 -0.02 -0.01 0.01 0.02 0.03 



t2 



Figure 1: Constraints on t2 and tiie linear combination (1 — j3) + from LLR and MRR, indicated 
by the arrows. For example, the hatched area is the region excluded by LLR only. General Relativity 
corresponds to /3 — 1 — t2 ^ ~ (black dot) . 



neglected, then simply > 0.0286 would be inconsistent with the data. Eventually, if 
(1 — /3 + 2t2 + ts) < —0.003 we would have analogous inconsistencies. 

We stress that the perihelion advance measurement used here is based on data taken between 
1966 and 1990. As pointed out by Will ([Tj, page 38) "analysis of data taken since 1990 
could improve the accuracy". Therefore, the above contraints on spacetime torsion can be 
improved already now, with existing data, while waiting that Mercury is reached by new 
spacecrafts, like in particular ESA's BepiColombo. 
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14 Equations of geodesies 



In this section we consider the particular case of geodesic trajectories. The system of 
equations of geodesies trajectories reads as 



dr"^ Ia^^J dr dr 



0. 



The resulting system of equations of motion is given by (6.3) with ^1 = ^2 = ^3 = 0. Such a 
system to lowest order becomes 

dv 1-Lm ^ 

Imposing the Newtonian limit (see also ^19^ formula (25)]) it follows that 



-2. 



Hence all the precession formulae are the same as in the PPN formalism given in (7.12), 



(11.5) and (11.9). Therefore, if a satellite's orbit is assumed to be a geodesic, then the 



measurements of satellite experiments cannot be used to constrain the torsion parameters. 



15 Discussion of the results and future prospects 

This work is an investigation of the effects of spacetime torsion on the orbits of satellites 
and planets, based on a model with several parameters evolved from the one by MTGC. 
Due to the presence of a set of parameters, it must be tested with a combination of ex- 
periments designed to measure different physical effects and observables. In this case, no 
single experiment provides a complete answer, but experiments with the best accuracy and 
the broadest parameter sensitivity may find the first reliable hint of torsion. The notable 
example is the geodetic precession, which can be measured using three very different in- 
strumental techniques: LLR, GPB's gyroscopes and future BepiColombo's radio science 
and accelerometer payloads. This makes constraining torsion with the geodetic precession 
robust against the effect of experimental systematic errors. 

For completeness, we quote here that the constraints on torsion provided by the measure- 
ment of the geodetic precession turn out to be useful also in constraining spacetime torsion 
with the frame dragging experiments on LAGEOS satellites. We remark that the torsion 
corrections to the Lense-Thirring effect for LAGEOS and GPB contain different sets of 
torsion parameters. We refer to the companion paper [12] for the details. 



15.1 Future prospects 

Before the end of the decade, robotic missions on the lunar surface could deploy new sci- 
entific payloads which include laser retroreflectors and thus extend the LLR reach for new 
physics in three ways: (i) using a significantly improved 2"'' generation retroreflector de- 
sign; (ii) increasing by a factor about 2 the geometric lever arm of LLR with missions to 
the lunar poles or limb; (iii) combining LLR payloads with transponders (at least two) for 
same-beam microwave interferometry (SBI) capable of additional accurate measurements 
of lunar rotations and librations [5^, [51] (during the lifetime of the transponders). 
In particular, the single, large, fused-silica retroreflector design developed by the University 
of Maryland and INFN-LNF [52j will improve over the performance of current Apollo arrays 
by a factor 100 or more, thus removing the dominant contribution to the LLR error budget. 



27 



Such a contribution is of the order of 2 cm. It is due to the multi-retroreflector structure 
of the arrays coupled to the hbrations and rotations of the Moon with respect to the 
Earth. The functionahty of this specific new design, which inherits and is evolved from 
the successful Apollo 11, 14 and 15 experience, is being validated by thermal- vacuum- 
optical testing laboratory-simulated space conditions at the INFN-LNF "Satellite/lunar 
laser ranging Characterization Facility (SCF)" [53], [5l], [55 j . 

Other instruments, like the seismometer and the heat-flow probe, will provide important 
information to evaluate the LLR systematic errors related to the environmental conditions 
of the lunar surface and sub-surface layers of the lunar regolith. 

After the end of this decade, results from the BepiColombo Mercury orbiter, an ESA Cor- 
nerstone mission equipped with a high-accuracy radio science and accelerometer payloads 
to test GR, is expected to improve the classical test of the perihelion advance f56], [57| 
(42.98 arcsec/century) . The latter measurement can be cross-checked by new MRR data 
taken simultaneously with BepiColombo's ranging data (possibily, by the same ground sta- 
tions). In addition, we note that in the past a Mercury orbiter like BepiColombo has 
been considered also for yet another independent measurement of the geodetic precession 
[58] , [59j (20.50 arcsec/century; to be compared to the 1.92 arcsec/century for the Moon). 
Mercury's special role in the search for new physics effects, and for spacetime torsion in 
particular, is due to the relatively large value of its eccentricity and to its short distance to 
the Sun. 

In conclusion, using current LLR and MRR data we have set constraints on the dimension- 
less torsion parameter t2 and the linear combination (1 — /?) + ^3 at 10^^ level. In the future, 
LLR, MRR, GPB and, ultimately, BepiColombo together can exclude non-zero values of t2 
and with accuracies significantly below 1%. 



16 Appendix 



In this appendix we compute the precession of pericenter using the method in Appendix C 
of the online version 1601. 



^/ \ dr dt 

+ (rV + rV.)--=o, 



The \ = t component of (6.1) is 

d^ 

which yields 



dH 
dr^ 

from which one finds 



where 



+ 



ti-n + {n'^ -21 + 2h) — 
exp I I il){r) dr 



A. 
d^ 



dt 
d^ 



m dr dt 
dr dr 

= 0, 



0, 



i^{r)= ti-n + {n^ -2X + 2t^) 
It follows the conservation law 



m 
r 



m 



k = exp 



1 



ti-n + -{n^ -21 + 2*3) 



ml dt 



r I dr 



const. 



The X = 4> component of (6.1) when = tt/2 is 

d\ 



dr^ 



dr dcp 
dr dr 



(16.1) 
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which yields 



(f(h f 2 m\ dr d(b 



^ / dr dr 



where we can neglect the term with the factor —2^4^. It follows the conservation law 



h = r exp t2 — 
dr \ r 



m 



const. 



(16.2) 



Notice that — 1 and h"^ are of the order em- Since the parameter r is the proper time we 
have 

ds^ dx^' dx" 



-1, 



dr'^ dr dr 

from which, for a test body in the equatorial plane 9 = ti/2 it follows 



, _ ,m 



dt 
dr 



m 



dr 



dr 



1. 



(16.3) 



Observe that the term with the factor I is missing in [60', (C23)]. 



Using (16.1), (16.2) and the identity dr/dr = {dr / d<j)){d<j)/ dr), from (16.3) it follows 



dr 



r4exp(2t2^) 



-1 + 



exp {-2 [t^-n + \{n^-2Z + 2U)f\ ^} 

h'' 1 



exp(2t2- 



(16.4) 



We now need to expand the right hand side of (16.4) up to the order e^, since the left hand 



side is of order 0(1). The right hand side is divided by h?, therefore we need to develop 
the quantity in {• • • } up to the order e^. Since /c^ — 1 is of the order e^, it is enough to 
develop the exponential in front of {• • • } up to the order em- Taking into account that h? 
is of order e^, we have 



/dry _ 



1 + (2t2 - F) 



m 



-l + k^ + k^ {2ti-U) 



m 



(16.5) 



[2ti{ti-n)+n'^ - i + 2t3] 



m 



h? om 

^ + 2*2^1% 



d2 



U 



From (16.5), setting u = 1/r, we obtain the differential equation of the orbit 
where 



+ u 



^ [k^ {-n -T + 2ti+ 2t2) +F- 2t2] + -Fmu^ + -^Am'^u, 



(16.6) 



(16.7) 



A = nT + W -T-2ti{n + T) + 2tl + 4tit2 - 2*2^ + 2*3, 

where we have neglected the term containing v?, which is of the order e^. 
We stress that the term j^Am'^u, neglected in the computations of [19], is of order Cm, as 
well as the term | J^mn^. Note also that ^ = in the case of GR, and that A = 2(2 — 7 — /?) 
in the case of PPN. 
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Using (6.5) and the above mentioned order of magnitude of k'^ it follows 



+ 1 



A 



I u 



m 



(16. 



In ( 16.8 ) we have neglected terms independent of u, originated by the expansion of k"^ inside 



in (16.6), which have no effect on the precession of the pericenter. 



If we neglect the last addendum on the right hand side of (16.8), and we set 

3 



m 



ci 



5i = -Tm, 



the solution uq of the corresponding linear equation is 



52 = j^m , 



(16.9) 



^o(0) = 1 |l + ecos yr 
1 - 02 ^ L 



ho)]} 



where e and 0o are constants of integration. If we neglect 82 with respect to 1, the above 
solution gives an elliptical orbit with eccentricity e and semi-latus rectum p = a(l — e^) = 



^— ^. Substituting Uq in place of in (16.8), the solution u of the corresponding linearized 



Cl 



equation is u = uq + ui, where 

5icl rl + eV2 



Ul{(f)) 



(1-52)' 



+ 



--{(t> - 0o) sin [yi - 82 (0 - </)o)) 



1 



l-Jo 6 



cos 



2^1 - 82 (0-00) 



If the eccentricity of the orbit is small (i.e., << e) and we neglect the (small) additive 



constant 
have 



u{(t)) ~ 



Sicj l+eV2 
(l-<52)2 l-<52 



which does not contribute to the precession of the pericenter, we 



Cl 



1 + e 



cos 



Taking into account that (-;^_^^^p/2 



n(0) 



Cl 



1-82 

Cl 



0) is of order and that 82 is of order em, we find 

, ^ 82 



1 



1 + e cos 
1 + e cos 



^0 - ci(5i(0 - ( 
0o) ( 1 - ci(5i 



(1 + 3ci<5i) (0 - 0o) 



8_2 

2 



This is the equation of an elliptic orbit whose pericenter precedes according to 



27r 



1 



1 - ci8 



1 - ^ 
1 2 



1 ~ 27r ci8 



(16.10) 



Substituting (16.7) and (16.9) into (16.10) we obtain formula (7.11). 
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